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We describe a peculiar fine structure acquired by the in-plane optical phonon at the T-point 
in graphene when it is brought into resonance with one of the inter-Landau-level transitions in 
this material. The effect is most pronounced when this lattice mode (associated with the G-band 
in graphene Raman spectrum) is in resonance with inter-Landau-level transitions =>• +, f and 
— , 1 0, at a magnetic field Bq ~ 30T. It can be used to measure the strength of the electron- 
phonon coupling directly, and its filling-factor dependence can be used experimentally to detect 
circularly polarized lattice vibrations. 

PACS numbers: 78.30.Na, 73.43.-f, 81.05.Uw 



In metals and semiconductors the spectra of phonons 
are renormalized by their interaction with electrons. 
Some of the best known examples include the Kohn 
anomaly [l[ in the phonon dispersion, which originates 
from the excitation/de-excitation of electrons across the 
Fermi level upon the propagation of a phonon through 
the bulk of a metal and a shift in the longitudinal opti- 
cal phonon frequency in heavily doped polar semiconduc- 
tors 0. However, despite the transparency of theoretical 
models the observation of such effects is often obscured 
by the difficulty to change the electron density in a mate- 
rial, whereas in semiconductor structures containing two- 
dimensional (2D) electrons the density of which can be 
varied, the influence of the latter on the phonon modes is 
weak due to a negligibly small volume fraction occupied 
by the electron gas. In this context, a unique opportunity 
arises in graphene-based field-effect transistors [3[ , where 
the density of carriers in an atomically thin film (mono- 
layer 0,(1,01 or a bilayer 0) can be continuously varied 
from 10 13 cm~ 2 p-type to 10 13 cm~ 2 n-type. Several Ra- 
man experiments have already been reported [1, Q where 
the variation of carrier density in graphene changes the 
optical phonon frequency, in agreement with theoretical 
expectations (lfl ITll. [l2| . 

When graphene is exposed to a quantizing magnetic 
field, its electronic spectrum quenches into discrete Lan- 
dau levels (LLs) [HI]. Then, the optical phonon energy in 
graphene may coincide with the energy of one of the inter- 
LL transitions , a condition known as magnetophonon 
resonance [3, Recently, Ando has suggested [l6| 
that in undoped graphene the magnetophonon resonance 
enhances the effect of the clcctron-phonon coupling on 
a spectrum of the in-plane optical phonons - the E 2s 
modes attributed to the G-band in the Raman spectra 
in Refs. [1, 0, [L7, 18, In this paper, we investigate 
a rich structure of the anti-crossing experienced by such 
lattice modes when a magnetic field makes their energy 
equal to the energy of one of the valley-antisymmetric 
interband magnetoexcitons [20| . Most saliently, the dif- 
ference between circular polarization of various inter-LL 



transitions 2l|, l22j makes the magnetophonon resonance 
distinguishable for lattice vibrations of different circular 
polarization, which makes the number of split lines in 
the fine structure acquired by a phonon and the value of 
splitting dependent on the electronic filling factor, v. 

The in-plane optical phonons in graphene [relative dis- 
placement u = (u x ,u y ) of sublattices A and B] have 
the energy ui « 0.2eV at the T-point (in the center of 
the Brillouin zone). These phonons and their coupling to 
electrons can be described using the Hamiltonian [lOllllI]. 



H Ph = y^^M.qVq + gVZMujVxUy - CTyUx), (1) 
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M , q are annihilation (creation) operators of a 
M is the mass of a car- 
is the number of unit cells. Here 
we use units ti = 1. Also, we shall uti- 



lize a double degeneracy of the E2 S mode at the T-point 
(at q = 0) and describe the in-plane optical phonon in 
terms of a degenerate pair of circularly polarized modes, 
uq = (u x + iu y )/^/2 and uq = u*q. The constant q in Eq. 

characterizes the electron-phonon coupling [23( . This 
coupling has the form of the only invariant linear in u per- 
mitted by the symmetry group of the honeycomb crystal. 
It is constructed using Pauli matrices er = (a x , o~ y ) acting 
in the space of sublattice components of the Bloch func- 
tions, [0k + a, 0k + b] and [</>k_b, </>k„a] which describe 
electron states in the valleys K± (two opposite corners of 
the hexagonal Brillouin zone) and obey the Hamiltonian, 
in terms of the electron charge — e < [2~j ], 



i?ei = S,vcr p, p = - iV + eA, d x A y 



d y A x = B. 



Here, £ = ± distinguishes between K±, and momentum 
p is calculated with respect to the center of the corre- 
sponding valley. This Hamiltonian represents the dom- 
inant term of the next-neighbor tight-binding model of 
graphene [HI, 0, [27| , and the electron-phonon coupling 
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FIG. 1: (a) Optical phonons are lattice vibrations with an cut- 
off-phase oscillation of the two sublattices. (b) Interband electron- 
hole excitations coupling to phonon modes with different circular 
polarization. 



in Eq. ([T]) takes into account the change in the A — B 
hopping elements due to the sublattice displacement 12811 . 

In a perpendicular magnetic field, H e \ determines |13f 
a spectrum of 4-fold (spin and valley) degenerate LLs, 



-a=± 
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in the valence band (£^ >0 )i con- 
duction band (£^>o)> an d at zero energy (eq = 0, ex- 
actly at the Dirac point in the electron spectrum), in 
terms of the magnetic length Xb = 1/veB. Such a 
spectrum has been confirmed by recent quantum Hall 
effect measurements 0, @, H|. In each of the two val- 
leys, the LL basis is given by two-component states 

\fkW l + <5n,o0n,m,i£a(l-<S n ,o)0n-i,m], where (j> n>m are 
the LL wave functions described by the quantum num- 
bers n and m, the latter being related to the guiding 
center degree of freedom. Here, we neglect the Zeeman 
effect, and simply take into account the two- fold spin de- 
generacy 

Excitations of electrons between LLs can be described 
in terms of magnetoexcitons (see Fig. 1). Those relevant 
for the magnetophonon resonance are 



E C +,n,m;t C -,{n+l),m;£, 



where the index A =0, characterizes the angular mo- 
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mentum of the excitation and the operators 
annihilate (create) an electron in the state a, n, m in 
the valley K^. The normalization factors jV„ = 
[(1 + 5 n , )NB(v-,( n +i) ~ v+.n)] 1 ' 2 and 7V„° = [(1 + 
8n,o)NB(i'-,n — v+^n+i))] 1 / 2 are used to ensure the 
bosonic commutation relations of the exciton operators, 

[>pA(n,0fi J A'( n '^')] = S A,A'h,e S n,n', where N B is the 
total number of states per LL in a sample, including 
the two-fold spin-degeneracy. These commutation rela- 
tions are obtained within the mean-field approximation 

witn ( c a.7um^ c a',n'.m':,e) = S^f6 a ,a'S n ,n' $m, m ' (S a - + 

5 a ,+^a.n), where < j> Qi „ < 1 is the partial filling fac- 



tor of the n-th LL. Similarly to magneto-optical selec 
tion rules in graphene 2^, 21. 22| 
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polarized phonons are coupled to electronic transitions 
with — , (n + 1) => +,n, and O-polarized phonons to 
—,n =>■ +, (n + 1) magneto-excitons, at the same en- 
ergy n„ = V2(7j/A B )(V^ + + 1) (Fig. [J), which 
follows directly from the composition of the LL in 
graphene and the form of the electron-phonon coupling 
in Eq. In contrast to photons that couple to 

the valley-symmetric mode ij) AjS (n) = [tp A (n, K+) + 
ip A (n, K_)]/\/2, electron-phonon interaction in Eq.JT]) 
couples phonons to the valley-antisymmetric magnetoex- 
citon ipA,as(n) = [ip A (n,K + ) - ip A (n,K-)]/\/2. 

In terms of magnetoexcitons we can, now, rewrite the 
electron-phonon Hamiltonian in a bosonized form, as 



H = E E^C^^rt^+E^ (3) 

r=s,as A.n A 

J2^( n ) b A tfj A;as (n) + b A ip A . as {n) , 
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(n+l) — v +,ni 



90 (n) = 9\J / (1 + Snfihy^ 



v +,{n+l), 



where g A are the effective coupling constants, with 7 = 
Z\/Za 2 /2irX B and a = 1.4A (distance between neighbor- 
ing carbon atoms). In the Hamiltonian ([3]), we have omit- 
ted electronic excitations with a higher angular momen- 
tum which do not couple to the in-plane optical phonon 
modes (e.g., n =>• n', with n' =/= nil). The dressed 
phonon operator corresponding to the Hamiltonian ([3]) 
is obtained by solving Dyson's equation. The pole of the 
propagator gives the antisymmetric coupled mode fre- 
quencies UJ Al 
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where n p stands for the number of the highest fully occu- 
pied LL in the spectrum, and A„ = >/2(v/\B){V n + 1 — 
\/n). In Eq. (|3|), the sum (extended up to the high- 
energy cut-off N ~ (Xb /a) 2 above which the electronic 
dispersion is no longer linear) takes into account inter- 
band magnetoexcitons, and the last term gives a small 
correction due to an intraband magnetoexciton. In the 
small- field limit and large doping (np 3> LL solution of 
Eq. reproduces the zero-field result [ljj, 03 if ° nc 
replaces the sum by an integral, Yln=o ~ > V dn, ap- 
proximates y/n. + \Jn + 1 ~ 2-y/n and A„ F f=a 0, and, 
then, linearizes Eq. Q by replacing Cu A by to in the de- 
nominator, 



uj ~ uiq + X 
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where A = (2/y/3Tr)(g/i) 2 ~ 3.3 x 10 3 is the same as in 
Refs. [13, Ei (t = 2v/3a ~ 3eV is the A — B hopping am- 
plitude), and (Do is the renormalized phonon frequency in 
an undopcd graphcne sheet at B = 0. The only variation 
arises at high fields, cDrj > ^/2v/Xb, where for np = the 
linearized Eq. (g]) yields 
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The strongest effect of the phonon coupling to elec- 
tron modes occurs when the frequency of the former 
coincides with the frequency Q n of one of the magne- 
toexcitons ^A.as{ n )- In such a case, the sum on the 
right-hand-side of the eigenvalue equation ^ is domi- 
nated by the resonance term and may be approximated 
by 2wg\(n)/ (u>_a ~ f^n)- This results in a fine structure 
of mixed phonon-magnetoexciton modes, ipA^sin) cos 9+ 
bj\ sin 9 with frequency Cu^ and tpA,as (n) sin 9 — cos 9 
with frequency (D^ [where cot 29 = (0„— u]o)/2g^\, which 
are determined for each polarisation (A =0, O ) sepa- 
rately, 



i (f2 n + fi)o) T y/^(n n -u> ) 2 +9^(n). (5) 



A generic form of the phonon-magnetoexciton anti- 
crossing and formation of coupled modes, aj^(n) in un- 
dopcd graphcne (i.e., v — 0) is illustrated in Fig. [2ja). 
Such an anticrossin g a nd mode mixing is simlar to that 
described by Ando [l(|. It can manifest itself in Raman 
spectroscopy: in a fine structure acquired by the G-line 
(earlier attributed [|, @, [H 0, 01] to the in-plane op- 
tical phonon at the T-point, E 2g mode) at the magneto- 
phonon resonance conditions. The effect is the strongest 
for the resonance £l n= o ~ ujq between the phonon and 
magnetoexciton based upon — , 1 => and =>• +, 1 tran- 
sitions. When approaching the resonance (by sweeping a 
magnetic field), the phonon line becomes accompanied by 
a weak satellite moving towards it and increasing its in- 
tensity. Exactly at the magnetophonon resonance, where 
both the upper mode [(I>j^(n)] and the lower mode [u)^(n)] 
consist of an equal-weight superposition of the phonon 
and the resonant exciton, with cos 6* = sin (9 = l/\/2, 
the G-band in graphene would appear as two lines. For 
n n=0 = V2v/X B ~ 36-/B[T] meV (see [H Q) and 
ujq ~ 200 meV, this resonance occurs in an experimen- 
tally accessible field range, Bo ~ 30 T. For the filling 
factor v = 0, the central LL (n = 0) is always half-filled. 
Then, coupling and, therefore, splitting of the O- and O- 
polarizcd modes coincide, g^ = go, thus, giving rise to a 
pair of peaks at the energies uj^ = loo ± go sketched 
in part I in Fig. [2][b). For the magnetic field value 
Bo — 30 T and g ~ 0.28eV [l2|, we estimate this splitting 
as 2gA ~ 16meV (~ 130cm _1 ), which largely exceeds the 
G-band width observed in Refs. @, 0, [H El El ■ 

Doping of graphene changes the strength of the cou- 
pling constants go and go, as shown in Fig. [^c). This 
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FIG. 2: (a) Coupled phonon and magneto-excitons as a function 
of the magnetic field. Energies are in units of the bare phonon 
energy w. Dashed lines indicate the uncoupled valley-symmetric 
modes, with g_4 = 0. (b) Mode splitting as a function of the filling 
factor, as may be seen in Raman spectroscopy, with the resonance 
condition Q n= o ~ wq, for u = in (I), < \v\ < 2 (in II), and 
v = ±2 (in III). The absolute intensity of the modes is in arbitrary 
units, but the height and the width reflect the expected relative 
intensities, (c) Mode splitting for n = 0, as a function of the filling 
factor v. (d) Same as in (c) for n > 1. 
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is because a higher (lower) occupancy of the n = LL 
reduces (enhances) the oscillator strength of the O polar- 
ized transition due to the availability of filled and empty 
states in the involved LLs. whereas the same change in 
the electron density has the opposite effect on g&. As 
a result, for an arbitrary filling factor —2 < v < 2, 
we predict that, in the vicinity of magnetophonon reso- 
nance, the phonon mode (and, therefore, G-band in Ra- 
man spectrum) should split into four lines [part II in Fig. 
EJb)], with u)q = tD±<7o for O-polarized and ojq = uj±go 
for O-polarized phonons. In the quantum Hall state at 
filling factor v = 2, the transition — , 1 =>■ becomes suc- 
cessively blocked and no longer affects the frequency of a 
O-polarized phonon, whereas the transition =>■ +, 1 ac- 
quires the maximum strength, thus, increasing the cou- 
pling parameter g^. This leads to the magnetophonon 
resonance fine structure consisting of three peaks, with 
an even larger splitting between side lines, as sketched in 
part III in Fig. EJb). Interestingly, this may enable one 
to directly observe lattice modes with a definite circu- 
lar polarization. A further increase of the electron filling 
factor reduces the side-line splitting which should com- 
pletely disappear at v = 6, after the transition =>■ +, 1 
becomes blocked by a complete filling of the +, 1 LL [Fig. 
[2fc)]. The same arguments hold for p-doped graphene, 
though in this case the roles of O- and O-polarized modes 
are interchanged. 

Magnetophonon resonances with other possible inter- 
LL transitions n n + 1 occur at much lower magnetic 
fields, B n = Bo/(^n+\/n + l) 2 . For example, a resonant 
phonon coupling with the magnetoexciton ipA;asO-) is ex- 
pected to occur at B\ f=a 5T. Its description remains qual- 
itatively similar, though for n > the mode splitting is 
less pronounced because of the _B-ficld dependence of the 
coupling constants in Eq. Q. One finds that g^j = g f j 
for \v\ < 2(2n - 1). At v = 2(2n - 1), filling of the n- 
th LL starts changing, which reduces splitting of the O- 
polarized mode and gives rise to the four-peak structure. 
At v = 2(2n+ 1), where the +, n LL becomes completely 
filled, splitting of the O-polarized phonon vanishes, thus, 
resulting in the three-peak fine structure [part III in Fig. 
EJb)] that would persist up to v = 2(2rt + 3). This is 
because the splitting of the O-polarized modes remains 
constant up to the filling factor v = 2{2n + 1), above 
which population of the +, (n + 1) LL starts to suppress 
the value of g&, until the latter vanishes at v = 2(2n + 3) 
[see Fig. m,d)]. 

In conclusion, we have predicted a filling-factor depen- 
dence of the fine structure acquired by the in-plane (E2 ff ) 
optical phonon in graphene when the latter is in reso- 
nance with one of the intcr-LL transitions in this ma- 
terial. The effect is expected to be most pronounced 
when the phonon is resonantly coupled to the +, 1 
and — , 1 transitions, which requires a magnetic field 
Bo ~ 30T. The predicted mode splitting may be used 
to measure directly the strength of the electron-phonon 



coupling, and also to distinguish between circularly (lcft- 
and right- hand) polarized lattice modes. 
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Erratum 

In the previous version (v3) of this Letter, we have 
underestimated the numerical value of the mode split- 
ting of the magnetophonon resonance [see paragraph af- 
ter Eq. (5)] by a factor of 2 (the text above takes into 
account the corrected parameters). This is a result of 
two mistakes. First, there is a factor of which finds 
its origin in an erroneous normalization of the circular 
polarized phonons. They should indeed be defined as 
Uq = (u x + iu y )j\f2 and u& = (u x — iu y )/\2 [and not 
as = u x + iu y and u& = u x — iu y as incorrectly 
assumed on page 1, second column], such that the asso- 
ciated phonon operators 6.4 obey the usual commutation 
relations [&a,6^,] = Sa,A', with A =0,0. This yields 
a factor of V2 in the definition of the effective coupling 
constants [Eq. (3)] , which read in the corrected form 



9o(n) = gy (1 + <5n,o)7-\/f-,(n+l) = v +,n , 
9o(n) = 9\J (1 + 8nfl)lyjv-,n ~ ^+,(n+l) ■ 

As a consequence, the zero-field dimensionless coupling 
constant A [defined in the first column page 3 of our 
Letter] is multiplied by a factor of 2 and becomes A = 
(2/V3n)(g/tf. 

Second, we also underestimated the numerical value of 
the electron-phonon coupling constant g by a factor of 
\[2. Indeed, g defined in our work [see Eq. (1)] is related 
to (9t)f — 0.0405 eV 2 computed by Piscanec et al. [|} 
as g = y / 2(g^) f — 0.28 eV and not as g = \J (g 2 ) p ~ 0.2 
cV as incorrectly assumed in our Letter. In addition, 
there is a substantial uncertainty in the precise value of 
the constant g. In a tight-binding model, the latter may 
be related to the derivative of the hopping amplitude t 
as a function of the carbon-carbon distance a as g — 
(—dt/da) x Harrison's phenomcnological 



law t oc 1/a 2 then implies that g ~ 0.26 cV. Experiments 
in graphene Q and Q in zero magnetic field give for the 
dimensionless coupling constant A the values 4.4 x ICO 3 
and 5.3 x 10 -3 respectively. This determines g in between 
0.3 eV and 0.36 eV, where we take into account that the 
value of t lies between 2.7 and 3 eV. In the end, we have 
to take g in the range between 0.26 and 0.36 eV [instead 
of g ~ 0.2 eV] and therefore the dimensionless coupling 
constant becomes A ~ (2.8 to 5.3) x 10~ 3 [instead of A ~ 
10" 3 ]. 

As a result of the two factors of the numerical 
estimate for the mode splitting 2gj^ at v = and B ~ 30 
T [at the discussed resonance — , 1 => and => +,1, see 
second column of page 3] becomes 2g_4 ~ 15 meV (~ 120 
cm" 1 ), for g ~ 0.26 eV and 2g A - 20 meV (~ 160 cm" 1 ) 
for g ~ 0.36 cV [instead of 2g A ~ 8 meV]. The effect 
is therefore twice larger than initially predicted. The 
conclusions of our work remain unaltered. 

We would like to thank C. Faugeras and M. Potemski 
for having drawn our attention on the underestimated 
value of the mode splitting. See also their recent preprint 
where they measure the magnetophonon resonance [5j. 
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